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This  study  presents  an  original  control  algorithm  for  a  hybrid  energy  system  with  a  renewable  energy 
source,  namely,  a  polymer  electrolyte  membrane  fuel  cell  (PEMFC)  and  a  photovoltaic  (PV)  array.  A  single 
storage  device,  i.e.,  a  supercapacitor  (ultracapacitor)  module,  is  in  the  proposed  structure.  The  main  weak 
point  of  fuel  cells  (FCs)  is  slow  dynamics  because  the  power  slope  is  limited  to  prevent  fuel  starvation 
problems,  improve  performance  and  increase  lifetime.  The  very  fast  power  response  and  high  specific 
power  of  a  supercapacitor  complements  the  slower  power  output  of  the  main  source  to  produce  the 
compatibility  and  performance  characteristics  needed  in  a  load.  The  energy  in  the  system  is  balanced  by 
d.c.-bus  energy  regulation  (or  indirect  voltage  regulation).  A  supercapacitor  module  functions  by  supply¬ 
ing  energy  to  regulate  the  d.c.-bus  energy.  The  fuel  cell,  as  a  slow  dynamic  source  in  this  system,  supplies 
energy  to  the  supercapacitor  module  in  order  to  keep  it  charged.  The  photovoltaic  array  assists  the  fuel  cell 
during  daytime.  To  verify  the  proposed  principle,  a  hardware  system  is  realized  with  analog  circuits  for 
the  fuel  cell,  solar  cell  and  supercapacitor  current  control  loops,  and  with  numerical  calculation  (dSPACE) 
for  the  energy  control  loops.  Experimental  results  with  small-scale  devices,  namely,  a  PEMFC  (1200  W, 
46  A)  manufactured  by  the  Ballard  Power  System  Company,  a  photovoltaic  array  (800  W,  31  A)  manu¬ 
factured  by  the  Ekarat  Solar  Company  and  a  supercapacitor  module  (100F,  32  V)  manufactured  by  the 
Maxwell  Technologies  Company,  illustrate  the  excellent  energy-management  scheme  during  load  cycles. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Currently,  most  of  the  energy  demand  in  the  world  is  met  by 
fossil  and  nuclear  power  plants.  A  small  part  is  drawn  from  renew¬ 
able  energy  technologies  such  as  wind,  solar,  fuel  cell,  biomass  and 
geothermal  energy  [1,2].  Wind  energy,  solar  energy  and  fuel  cells 
have  experienced  a  remarkably  rapid  growth  in  the  past  ten  years 
[3-5]  because  they  are  pollution-free  sources  of  power.  Addition¬ 
ally,  they  generate  power  near  the  load  centres,  which  eliminates 
the  need  to  run  high-voltage  transmission  lines  through  rural  and 
urban  landscapes. 

The  cost  of  solar  photovoltaic  and  fuel  cell  electricity  is  still  high 
[6-8].  Nevertheless,  with  ongoing  research,  development  and  uti¬ 
lization  of  these  technologies  around  the  world,  the  costs  of  solar 
cells  and  fuel  cell  energy  are  expected  to  fall  in  the  next  few  years.  As 
for  solar  cell  and  fuel  cell  electricity  producers,  they  now  sell  power 
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freely  to  end-users  through  truly  open  access  to  the  transmission 
lines.  For  this  reason,  they  are  likely  to  benefit  as  much  as  other 
producers  of  electricity.  Another  benefit  in  their  favour  is  that  the 
cost  of  renewable  energy  falls  as  technology  advances,  whereas  the 
cost  of  electricity  from  conventional  power  plants  rises  with  infla¬ 
tion.  The  difference  in  their  trends  indicates  that  hydrogen  and  solar 
power  will  be  more  advantageous  in  future. 

In  the  near  future,  the  utility  power  system  at  a  large  scale  will 
be  supplied  by  renewable  energy  sources  and  storage  device(s),  i.e., 
hybrid  energy  systems ,  in  order  to  increase  their  reliability  and  make 
them  more  effective.  The  specific  properties  of  fuel  cells  and  solar 
cells  are  as  follows: 


•  The  amount  of  power  the  fuel  cell  system  delivers  is  controlled 
by  the  amount  of  current  drawn  from  the  system  if  the  proper 
conditions  for  cell  operation  are  maintained.  When  a  large  load 
is  applied  to  the  cells,  the  sudden  increase  in  the  current  can 
cause  the  system  to  stall  if  the  depleted  oxygen  or  hydrogen 
cannot  be  replenished  immediately  and  sufficiently.  Cell  star¬ 
vation  can  lead  to  a  system  stall,  permanent  cell  damage  or 
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Nomenclature 

a.c. 

alternating  current 

d.c. 

direct  current 

FC 

fuel  cell 

PV 

photovoltaic 

SC 

supercapacitor 

Cbus 

total  capacitance  at  d.c.-bus  (F) 

Csc 

total  capacitance  of  supercapacitor  module  (F) 

*Load 

d.c.-bus  load  current  (A) 

i'pc 

fuel  cell  current  (A) 

*FCREF 

fuel  cell  current  reference  (set-point)  (A) 

*Sol 

solar  cell  (photovoltaic)  current  (A) 

isolREF 

solar  cell  current  reference  (set-point)  (A) 

he 

supercapacitor  current  (A) 

*SCREF 

supercapacitor  current  reference  (set-point)  (A) 

PLoad 

load  power  (W) 

Pfc 

fuel  cell  power  (W) 

PFCa 

fuel  cell  output  power  to  d.c.-bus  (W) 

Pfcref 

fuel  cell  power  reference  (set-point)  (W) 

PFCMax 

maximum  fuel  cell  power  (W) 

PSol 

solar  cell  (photovoltaic)  power  (W) 

PSola 

solar  cell  output  power  to  d.c.-bus  (W) 

PSolREF 

solar  cell  power  reference  (set-point)  (W) 

PSolMax 

maximum  solar  cell  power  (W) 

Psc 

supercapacitor  power  (W) 

PsCa 

supercapacitor  output  power  to  d.c.-bus  (W) 

PSCREF 

supercapacitor  power  reference  (set-point)  (W) 

PsCMax 

maximum  supercapacitor  power  (W) 

^Bus 

d.c.-bus  voltage  (V) 

I'PC 

fuel  cell  voltage  (V) 

^Sol 

solar  cell  (photovoltaic)  voltage  (V) 

VSC 

supercapacitor  voltage  (V) 

TBus 

d.c.-bus  energy  (J) 

yBusREF 

d.c.-bus  energy  reference  (set-point)  (J) 

ysc 

supercapacitor  energy  (J) 

ySCREF 

supercapacitor  energy  reference  (set-point)  (J) 

yr 

total  energy  at  d.c.-bus  and  supercapacitor  (J) 

rFc 

equivalent  series  resistance  in  fuel  cell  converter 

m 

rs  oi 

equivalent  series  resistance  in  solar  cell  converter 

m 

rsc 

equivalent  series  resistance  in  supercapacitor  con¬ 
verter  (£2) 

u 

input  variable  vector 

X 

state-variable  vector 

y 

output  vector 

<p(-x  n- 

),  0(  • )  smooth  mapping  functions 

reduced  cell  lifetime.  To  protect  the  fuel  cells  from  overloading 
and  starvation,  especially  during  transient  conditions,  excessive 
oxygen  and  hydrogen  can  be  supplied  to  the  cells  during  the 
steady-state  operation,  which  increases  the  reserve  of  available 
power  in  anticipation  of  a  load  increase.  This  strategy,  however, 
is  conservative  and  leads  to  increased  parasitic  losses,  decreased 
air  utilization  and  thereby  compromised  system  performance. 
Therefore,  the  fuel  cell  power  or  current  slope  must  be  limited 
to  prevent  a  fuel  cell  stack  from  experiencing  the  fuel  starvation 
phenomenon  and  to  optimize  the  system,  for  example,  4  A  s_1  for 
a  0.5 kW,  12.5-V  PEMFC  [9],  and  5,  10  and  50 As-1  for  a  20 kW, 
48 -V  PEMFC  [10]. 

•  The  output  power  of  solar  cells  fluctuates  considerably  depending 
on  solar  radiation,  weather  conditions  and  temperature  [11,12]. 


Therefore,  in  order  to  supply  electric  power  to  fluctuating  loads 
with  a  hybrid  system  composed  of  the  above-mentioned  fuel  cell 
and  solar  cell,  an  electric  energy-storage  system  is  needed  to  com¬ 
pensate  for  the  gap  between  the  output  from  the  renewable  energy 
sources  and  the  load,  in  addition  to  the  collaborative  load  sharing 
among  those  energies  [13-15]. 

Moreover,  hydrogen  as  an  energy  storage  media  has  the  poten¬ 
tial  to  address  both  daily  and  seasonal  buffering  requirements. 
Systems  that  employ  an  electrolyzer  to  convert  excess  electric¬ 
ity  to  hydrogen  coupled  with  hydrogen  storage  and  regeneration 
using  a  fuel  cell  can,  in  principle,  provide  power  with  zero  (or 
near  zero)  emissions.  Hydrogen  production  by  solar  energy  is  a 
‘renewable-regenerative  system’  [4],  and  this  process  is  known  as 
the  electrolysis  process.  The  basic  principle  is  the  following:  when 
the  photovoltaic  input  power  exceeds  the  load  power  demand,  the 
system  controller  determines  that  the  energy  should  be  directed 
to  hydrogen  production.  In  this  kind  of  operation,  i.e.,  a  solar- 
based  renewable-regenerative  system,  almost  half  of  the  solar 
input  energy  is  directed  to  hydrogen  production  and  converted 
with  60%  energy  efficiency  [4]. 

Based  on  present  storage  device  technology,  battery  design  has 
to  supply  the  trade-off  between  specific  energy,  specific  power  and 
cycle  life.  The  difficulty  in  obtaining  high  values  of  these  three 
parameters  has  led  to  some  suggestions  that  the  energy-storage 
system  of  distributed  generation  systems  should  be  a  hybridiza¬ 
tion  of  an  energy  source  and  a  power  source  [16,17].  The  energy 
source,  mainly  fuel  cells  and  solar  cells  in  this  study,  has  high  spe¬ 
cific  energy,  whereas  the  power  source  has  high  specific  power.  The 
power  sources  can  be  recharged  from  the  main  energy  source(s) 
when  there  is  less  demand.  The  power  source  that  has  received 
wide  attention  is  the  supercapacitor  (or  ‘ultracapacitor’,  or  ‘elec¬ 
trochemical  double-layer  capacitor’)  [18-20]. 

The  enhancements  in  the  performance  of  renewable  energy 
source  power  systems  that  are  gained  by  adding  energy  storage 
are  all  derived  from  the  ability  to  shift  the  system  output.  Firming- 
up  the  renewable  system  is  accomplished  by  ensuring  that  energy 
is  available  when  there  is  a  demand  for  it  rather  than  being  limited 
by  the  availability  of  the  renewable  resource.  As  a  result,  the  sys¬ 
tem  output  may  need  to  be  shifted  to  periods  when  the  hydrogen 
and/or  the  sun,  for  example,  are  not  available  [21].  Depending  on 
the  size  and  type  of  the  energy-storage  system  and  the  load,  it  may 
be  possible  to  provide  all  of  the  power  needed  to  support  the  load. 
A  much  more  common  scenario  is  for  the  energy  storage  to  simply 
provide  enough  power  for  applications,  like  peak  shaving,  without 
having  full  load  support  capability.  The  energy-storage  system  may 
also  provide  sufficient  energy  to  ride  out  electric  service  interrup¬ 
tions  that  range  from  a  few  seconds  to  a  few  hours.  This  is  especially 
important  for  service  disruptions  that  occur  when  the  renewable 
resource  is  not  available. 

In  this  study,  a  fuel  cell/photovoltaic/supercapacitor  hybrid 
power  source  is  proposed.  A  power  electronic  converter  struc¬ 
ture,  energy  management  and  innovative  energy  control  law  are 
presented.  To  authenticate  the  proposed  principle,  the  hardware 
system  is  realized  by  analog  circuits  and  digital  estimation  with  a 
dSPACE  controller.  Experimental  results  obtained  with  small-scale 
devices  illustrate  the  system  performance. 


2.  Renewable  energy  hybrid  system 

2.1.  System  configuration 

The  power  converter  structure  of  the  system  is  shown  in  Fig.  1 . 
Fuel  cell  and  solar  cell  power-generating  systems  may  provide 
direct  or  alternating  current  (d.c.  or  a.c.)  to  satisfy  application- 
specific  power  needs  (Fig.  1 ).  The  current,  voltage  and  power  quality 
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Solar  Cell 

Solar  Cell  Panels  Converter 


Supercapacitor 


Fig.  1.  Proposed  hybrid  energy  system  supplied  by  fuel  cell,  solar  cell  and  supercapacitor,  where  pLoad  (=  vbus  x  iLoad).  vbus  and  iLoad  are  load  power,  d.c.-bus  voltage  and 
d.c.-bus  load  current,  respectively.  pFC  (=  vFC  x  Ifc),  vfc  and  iFC  are  fuel  cell  power,  voltage  and  current,  respectively.  pSo i  (=  Vsoi  x  *soi)>  vSoi  and  iSo i  are  solar  cell  power,  voltage 
and  current,  respectively.  pSc  (=  Use  x  isc),  Vsc  and  i sc  are  supercapacitor  power,  voltage,  and  current,  respectively.  pFca,  Psoia  and  pSc a  are  the  output  powers  to  d.c.  link 
from  converters  of  fuel  cell,  solar  cell  and  supercapacitor,  respectively.  i'FCref.  i’soiref.  i’scref  are  the  current  references  of  the  fuel  cell,  solar  cell  and  supercapacitor  converters, 
respectively. 


are  controlled  by  electronic  power-conditioning  systems.  Gener¬ 
ally,  voltage  regulators  and  d.c.-d.c.  converters  (chopper  circuits) 
are  used  to  control  and  adjust  the  fuel  cell  and  photovoltaic  d.c. 
output  voltages  to  useful  values  [22,23]. 

The  voltages  of  the  fuel  cell  and  solar  cell  vary  with  the  cur¬ 
rent  drawn  by  the  load  and  decrease  significantly  at  high  load 
currents  [24].  Therefore,  a  power-electronics  system  is  needed  to 
process  the  raw  output  power  from  the  stack  or  panel  and  pro¬ 
vide  power  to  the  load  at  a  constant  d.c.  or  a.c.  voltage.  Typically, 
the  power-electronics  systems  include  multiple  interconnected 
power  converters:  for  example,  a  d.c.-d.c.  converter  is  often  fol¬ 
lowed  by  a  d.c.-a.c.  converter  for  stationary  applications  [25].  The 
switching  scheme  in  such  power  converters  can  be  based  on  pulse- 
width  modulation  (PWM),  resonant,  quasi-resonant,  soft-switched, 
or  line-commutated  switching  [26-28].  Furthermore,  the  topolog¬ 
ical  structures  of  these  converters  can  vary  dramatically  [29]. 

The  supercapacitor  bank  is  always  connected  to  the  d.c.-bus 
by  means  of  a  two-quadrant  d.c.-d.c.  converter  (bidirectional  con¬ 
verter),  a  ‘supercapacitor  converter’.  The  supercapacitor  power  can 
be  positive  or  negative,  which  allows  energy  to  be  transferred  in 
both  directions  [30]. 

An  inverter  is  used  to  convert  the  d.c.  voltage  to  a  useful 
a.c.  voltage.  It  may  operate  as  grid-independent  or  with  grid- 
synchronization.  Then,  the  solar  cell  supplies  all  its  output  power 
to  the  d.c.-bus,  and  the  fuel  cell  supplies  a  slowly  fluctuating  power 
difference  between  the  load  power  pLoad  demand  and  the  solar  cell 
output  power.  A  supercapacitor  storage  device  can  compensate  for 
both  a  sudden  change  in  the  load  demand  and  rapidly  fluctuating 
solar  cell  output  power  caused  by  solar  radiation;  therefore,  the 
load  is  supplied  with  stable  power. 

As  depicted  in  Fig.  1,  the  proposed  system  has  a  multi-source 
structure  (or  cascade  scheme).  The  cascaded  structure  can  lead  to 
interactions  between  converters  that  cause  the  system  to  be  unsta¬ 
ble,  if  they  are  designed  separately.  This  problem  has  been  studied 
during  the  past  decade.  In  this  case,  to  study  interactions  between 
the  converters,  impedance  criteria  are  often  used  to  determine  the 


stability  of  the  cascaded  system.  Nonetheless,  the  technique  only 
proves  the  asymptotic  stability  around  a  given  operating  point.  The 
large  signal  stability  properties  or  the  behaviour  of  the  system  in 
the  case  of  outsized  disturbances  are  not  considered  in  this  model. 
In  this  paper,  a  non-linear  control  algorithm  based  on  the  differen¬ 
tial  flatness  principle  of  the  renewable  power  plant  is  proposed. 
Design  controller  parameters  are  independent  of  the  operating 
point.  Interactions  between  the  converters  are  taken  into  account 
by  the  controllers,  and  high  dynamics  in  perturbation  rejection  are 
accomplished. 

2.2.  Energy  management 

For  reasons  of  safety  and  dynamics,  the  fuel  cell,  solar  cell  and 
supercapacitor  converters  are  usually  controlled  primarily  by  inner 
current- regulation  loops  [17].  The  dynamics  of  current-regulation 
loops  are  much  faster  than  those  of  outer  control  loops,  which  are 
described  later.  Therefore,  the  fuel  cell  current  ipc,  the  solar  cell 
current  iSol  and  the  supercapacitor  current  iSc  are  considered  to 
follow  their  references  i'fcref.  i’soiref  and  i'scref  perfectly. 

Energy  management  of  multi-power  sources  has  already  been 
studied  recently,  for  example,  Feroldi  et  al.  [31]  studied  control 
(based  on  efficiency  map)  of  a  fuel  cell/supercapacitor  hybrid  power 
source  for  vehicle  applications,  Payman  et  al.  [32]  studied  control 
of  a  regulated  d.c.  voltage  supplied  by  a  fuel  cell  and  superca¬ 
pacitor  storage  device  based  on  a  differential  flatness  system, 
and  Thounthong  et  al.  addressed  a  regulated  d.c.-bus  voltage  fuel 
cell/supercapacitor  hybrid  source  [33],  a  regulated  d.c.-bus  voltage 
fuel  cell/battery/supercapacitor  hybrid  source  [34]  and  an  unregu¬ 
lated  d.c.-bus  voltage  fuel  cell/battery  hybrid  source  [14]. 

Thus,  in  the  proposed  system  depicted  in  Fig.  1,  there  are  two- 
voltage  variables,  or  two-energy  variables,  to  be  regulated.  The 
d.c.-bus  energy  yBus  is  the  most  important  variable,  and  the  super¬ 
capacitor  storage  energy  ySc  is  the  next  most  important.  Therefore, 
based  on  the  previous  work  described  above,  it  is  proposed  to  utilize 
supercapacitors,  which  are  the  fastest  energy  source  in  the  pro- 
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posed  system,  to  supply  the  energy  for  the  d.c.-bus  [33,34].  Thus, 
the  fuel  cell  (as  the  slowest  dynamic  device)  and  solar  cell  supply 
energy  for  both  the  d.c.-bus  capacitor  CBus  and  the  supercapacitors 
Csc  to  keep  them  charged. 

At  the  d.c.  link,  there  are  three  sources  that  supply  energy  at  the 
same  time.  Therefore,  the  dynamics  of  the  three  sources  must  be 
classified  in  order  to  guarantee  system  stability  [34],  as  illustrated 
in  Fig.  2.  Note  that  the  power  unit  in  Fig.  2  is  the  normalized  unit  (or 
per-unit  [pu]).  Thus,  the  defined  dynamic  classification  depicted  in 
Fig.  2  is  obtained  from  the  simulation  result  using  Matlab-Simulink, 
and  the  experimental  results  of  the  power  source  dynamics  are  pre¬ 
sented  after  the  performance  validation  section.  The  supercapacitor 
is  operated  as  the  highest  dynamic  power  source  to  provide  the 
micro-cycles  and  the  fast  dynamic  power  supply.  It  can  withstand 
a  very  large  number  of  charge-discharge  cycles  without  degrada¬ 
tion  (or  virtually  infinite  cycles)  [33,35].  The  fuel  cell  generator  is 
operated  as  the  lowest  dynamic  power  source.  The  fuel  cell  cur¬ 
rent  or  power  slope  must  be  limited  to  avoid  the  fuel  starvation 
phenomenon.  The  fuel  cell  limited  current  or  power  slope  has  been 
experimentally  determined  to  be  the  highest  slope  of  an  operating 
fuel  cell  system,  where  no  fuel  starvation  occurs  in  order  to  improve 
its  lifetime  [33,34].  The  photovoltaic  generator  is  between  the  fuel 
cell  and  supercapacitor  in  the  dynamic  classification. 

2.2.1.  Principle  of  differential  flatness  theory 

A  non-linear  control  algorithm  based  on  the  flatness  properties 
[36-38  ]  of  the  system  is  proposed.  According  to  the  flatness-control 
law,  the  tuning  controller  parameters  are  independent  of  the  oper¬ 
ating  point  [32].  According  to  Fliess  et  al.  [39,40],  an  independent 
dynamic  system  is  considered  to  be  differentially  flat  if  an  output 
y  is  given  by: 

y  =  0(x,  u,  u,  ...,u(a)),  yeffim  (1) 

such  that  the  state  x  and  controls  u  (or  inputs)  can  be  written  as: 

x  =  <p(y,y,  ...,y(^),  (2) 

u  =  i//(y,y,  ...,y(^+1)),  ue9tm  (3) 

where  <p(  ),  0(-)  are  the  smooth  mapping  functions,  a  and  (3  are 

a  finite  number  of  their  time  derivatives,  and  (n,  m )  e  N.  Moreover, 
it  is  assumed  that  m  <  n. 

The  output  y  is  called  a  flat  output.  Thus,  a  dynamic  system  is, 
in  nature,  differentially  flat  if  it  is  equivalent  to  a  system  without 
dynamics,  i.e.,  a  static  system.  In  other  words,  there  are  no  dif¬ 
ferential  constraints  in  output  space.  The  apparent  advantage  of  a 


differentially  flat  system  is  that  every  trajectory  in  the  output  space 
is  feasible;  therefore,  trajectory  generation  is  theoretically  simpler 
in  terms  of  the  flat  outputs.  For  this  reason,  by  using  the  flatness 
estimation  in  power  electronic  applications,  interactions  between 
the  converters  are  taken  into  account  by  the  flatness-control  law 
and  high  dynamics  in  perturbation  elimination  are  achieved. 


2.2.2.  d.c.-Bus  energy  control 

Based  on  the  flatness-control,  the  electrostatic  energy  yBus 
stored  in  the  total  d.c.-bus  capacitor  CBus  is  considered  as  the  flat 
output  variable  (Eq.  (1)),  i.e.: 


1  7 

yBus  —  2^Bus1;Bus  (4) 

It  is  assumed  that  the  supercapacitor  current  follows  its  reference 
value  perfectly.  Thus: 


isc  =  *scref  =  — - 
vsc 


P SCREF 
VSC 


(5) 


where  Pscref  is  the  control  variable  of  the  reduced  system  described 
in  Eq.  (3). 

Therefore,  the  state  variable  i;Bus  may  be  put  into  a  form  similar 
to  that  in  Eq.  (2): 


^Bus 


=  <p\iy) 


(6) 


To  obtain  a  relationship  between  the  control  variable  pSc  and  the 
flat  output  variable  yBus,  the  differential  equation  verified  by  the 
reduced  model  from  Fig.  1  is  written  as: 


yBus  =  PFCa  +  PSola  +  PsCa  ~  PLoad 


(7) 


where 

PFCa=PFC-rFC(^) 
PSola  =  PSol  -  rsol  (^) 
Psca=Psc-rsc(g) 


(8) 

0) 

(10) 


rFC  is  the  total  equivalent  series  resistance  in  the  fuel  cell  converter. 
Here,  only  static  resistance  loss  is  considered  [32].  rSo \  is  the  total 
equivalent  series  resistance  in  the  solar  cell  converter.  rs c  is  the 
total  equivalent  series  resistance  in  the  supercapacitor  converter. 
It  follows  that: 


Pscref  =  2pScMax 

\  J  ( ^Bus  +  PLoad  -  PFCa  “  PSola  \ 

V  V  PsCMax  J 

=  f\(y,y)  (ii) 


where 


PsCMax  = 


y2 

VSC 

4rsc 


(12) 


PLoad  —  ^Bus^Load  — 


2yBu 


Cbu 


hoad 


(13) 


PscMax  is  the  limited  maximum  power  from  the  supercapacitor 
module. 

Thus,  it  is  apparent  thatpSc  =  iAi(y,  y),  which  corresponds  with 
Eq.  (3),  and  the  proposed  reduced  system  with  the  control  variable 
Psc  can  be  considered  to  be  a  flat  system  associated  withyBus  as  the 
flat  output  variable. 

To  ensure  the  control  of  the  flat  output  variable  yBus  to  its  ref¬ 
erence  trajectory  yBusREF,  the  following  control  law  is  based  on  the 
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Fig.  3.  Proposed  d.c.-bus  energy  control  loop  for  fuel  cell/solar  cell/supercapacitor  hybrid  source. 


well-known  second-order  control  law: 


(J/Bus  -  yBusREF)  +  ^11  (yBus  -  yBusREF) 


+f<12  /  (yB  us  —  yBusREF)  —  0 


(14) 


where  J<n,  and  K12  are  the  controller  parameters.  Defining 
eyBus  =yBus  -yBusREF,  I<n  =  2£ o)n ,  and  /Ci2  =  gives: 


^yBus  +  2^CdnCyBus  +  ^n^yBus  —  0 


(15) 


It  is  obvious  that  the  control  system  is  stable  for  I<u,  JC12>  0. 
Based  on  the  power  electronic  constant  switching  frequency /s  and 
a  cascade  control  structure,  the  outer  control  loop  (here  the  d.c.-bus 
energy  control)  must  operate  at  a  cut-off  frequency /e  <<c/c  (a  cut¬ 
off  frequency  of  the  supercapacitor  current  loop  or  power  loop)  </s 
[41  ].  The  d.c.-bus  energy  control  loop  detailed  above  is  portrayed 
in  Fig.  3. 

The  d.c.-bus  energy  control  law  generates  a  supercapacitor 
power  reference  Pscref-  This  signal  is  then  divided  by  the  measured 
supercapacitor  voltage  v$c  and  limited  to  maintain  the  superca¬ 
pacitor  voltage  within  an  interval  [VSCMin,  VscMax]  by  limiting  the 
supercapacitor  charging  current  or  discharging  current,  as  pre¬ 
sented  in  the  block  “SuperC  Current  Limitation  Function”  [33].  The 
higher  voltage  value  VSCMax  of  this  interval  corresponds  to  the  max¬ 
imum  voltage  of  the  storage  device.  Generally,  the  lower  voltage 
value  VSCMin  is  chosen  as  VscMax/2,  which  causes  only  25%  of  the 
energy  in  the  supercapacitor  bank  to  remain  [33];  as  a  result,  the 
supercapacitor  discharge  becomes  ineffective.  This  results  in  the 
supercapacitor  current  reference  i’scref- 


2.2.3.  Supercapacitor  energy  control 

The  classical  electrostatic  energy  storage  ysc  in  the  supercapac¬ 
itor  is  expressed  by: 

ysc  =  2Csci;sc  (16) 

where  CSc  is  the  supercapacitor  capacitance  (refer  to  Fig.  1 ),  and  vSc 
is  the  supercapacitor  voltage. 

Again,  based  on  the  flatness-control,  the  total  electrostatic 
energy  yT  stored  in  the  total  d.c.-bus  capacitor  CBus  and  the  super¬ 
capacitor  CSc  is  taken  as  the  flat  output  variable  (Eq.  (1 )),  i.e.: 

1  ?  1  ? 

yT  =  yBus  +ysc  =  2  CBusyBus  +  2  ^SC^SC  (I7) 


It  is  again  assumed  that  the  fuel  cell  and  photovoltaic  currents 
follow  their  reference  values  perfectly.  Thus: 


Pfc 

1  fc  =  *fcref  =  — - 
VFC 


lSol  =  *SolREF  : 


PSol 

I'Sol 


Pfcref 

VFc 

_  PSolREF 
I'Sol 


(18) 

(19) 


where  Pfcref  and  PsoIref  are  the  control  variables  of  the  reduced 
system  (Eq.  (3)). 


Thus,  the  state  variable  vsc  may  be  put  into  a  form  similar  to  that 
of  Eq.  (2): 

Vsc  =  ]p^^  =  (p2{y)  (20) 

The  differential  equation  verified  by  the  reduced  model  from  Fig.  1 
is: 

yT  =  PTa-PLoad  (21) 

where 

PTa  =  PFCa  +  Psola  =  PT  —  *T  (  ~ )  (22) 

For  a  fuel  cell  and  a  solar  cell  that  are  combined  as  an  energy  source 
to  supply  energy  to  a  d.c.-bus  and  supercapacitor,  vj  is  defined  as 
presented  in  Eq.  (22)  and  rT  represents  the  total  losses  in  the  fuel 
cell  and  photovoltaic  converters. 


Therefore: 

Ptref  =  2pTMax 

'-f- 

/Tt  +  PLoad  \ 

V  PTMax  J 

(23) 

where 

V2J 

p™ax  =  4^ 

(24) 

In  fact: 

PTMax  =  PFCMax  +  PFCMax 

(25) 

PFCMax  is  the  maximum  fuel  cell  power.  PsoiMax  is  the  maximum 
photovoltaic  power. 

Photovoltaic  power  systems  require  some  specific  estimation 
algorithms  to  deliver  the  maximum  power  point  (MPP)  because 
the  output  features  of  a  solar  cell  vary  with  environmental  changes 
in  irradiance  and  temperature.  Many  previous  investigations  have 
presented  different  techniques  for  maximum  power  point  tracking 
(MPPT)  of  photovoltaic  arrays  [42-44].  It  is  beyond  the  scope  of  this 
study  to  present  the  MPPT  algorithm  here,  though  it  is  similar  to  Eq. 
( 14).  Because  the  supercapacitor  energy  storage  has  a  large  capacity 
and  the  supercapacitor  energy  is  defined  as  a  slower  dynamic  vari¬ 
able  than  the  d.c.-bus  energy  variable,  the  total  (supercapacitor) 
energy  control  law  is  defined  as: 

(yT  -  J/TREf)  +  J<21  (yT  -  yTREF )  =  0  (26) 

The  supercapacitor  energy  control  loop  is  shown  in  Fig.  4.  The 
total  energy  control  law  generates  a  total  power  reference  Ptref- 
First,  Ptref  is  considered  to  be  the  photovoltaic  power  reference 
PsoiarREF-  It  must  be  limited  in  level,  within  an  interval  maximum 
PsoiarMax  (MPP)  and  minimum  pSoiarMin  (set  to  0  W),  and  have  lim¬ 
ited  dynamics  to  respect  constraints  in  the  dynamic  classification 
(Fig.  2).  Here,  the  second-order  delay  [34]  is  selected  as  a  photo¬ 
voltaic  power  dynamic  limitation,  as  illustrated  in  Fig.  4.  Second, 
the  difference  between  the  total  power  reference  Ptref  and  the 
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Fig.  4.  Proposed  supercapacitor  energy  control  loop  for  fuel  cell/solar  cell/supercapacitor  hybrid  source. 


photovoltaic  power  reference  PsoiarREF  is  the  fuel  cell  power  refer¬ 
ence  Pfcref-  It  must  be  limited  in  level,  within  an  interval  maximum 
PFCMax  and  minimum  pFCMin  (set  to  0  W),  and  be  limited  in  dynam¬ 
ics  to  respect  the  dynamic  classification  constraints  in  Fig.  2.  Again, 
second-order  delay  is  selected  as  a  fuel  cell  power  dynamic  limita¬ 
tion. 

3.  Performance  validation 

3.1.  Test  bench  description 

In  order  to  authenticate  the  proposed  original  hybrid  energy 
management,  a  small-scale  test  bench  has  been  implemented,  as 
presented  in  Fig.  5.  The  fuel  cell  system  used  in  this  effort  is  a  Bal¬ 
lard  Nexa  fuel  cell  system  (1.2  kW,  46  A).  It  is  supplied  using  pure 
hydrogen  with  a  regulated  pressure  of  lObar  from  bottles  under 
a  pressure  of  150  bar  and  with  clean  dry  air  from  a  compressor. 
The  photovoltaic  array  (800  W,  31  A)  was  obtained  by  connecting 
four  panels  (200  W,  7.78  A)  in  series.  The  solar  cell  panels  were 
developed  and  manufactured  by  the  Ekarat  Solar  Company.  They 
were  installed  on  the  roof  of  the  laboratory  building,  as  illus¬ 
trated  in  Fig.  5.  The  supercapacitor  module  (100F,  32  V,  based  on 
Maxwell  Technologies  Company)  was  obtained  by  connecting  12 
cells  BCAP1200  (capacitance:  1200  F;  maximum  voltage:  2.7  V)  in 
series,  as  shown  in  Fig.  5. 

The  fuel  cell  converter  (1200W)  is  composed  of  four-phase 
boost  converters  connected  in  parallel  using  the  interleaving  tech¬ 
nique  [45]  and  is  selected  to  adapt  the  low  d.c.  voltage  delivered 


by  the  fuel  cell,  which  is  approximately  26  V  at  rated  power,  to  the 
60  V  d.c.-bus.  Thus,  each  boost  converter  is  composed  of  a  high- 
frequency  input  inductor  (420  fxH),  four  output  filtering  capacitors 
connected  in  parallel  (each  capacitor  is  an  aluminum  electrolytic 
capacitor  (470  |xF),  thus  the  total  capacitance  is  1880  pT),  a  diode 
(Schottky  diode  RURG3020:  200V,  30 A),  and  a  main  power 
switch  (power  metal-oxide-semiconductor  field-effect  transistor 
IRFP264N  MOSFET:  250  V,  38  A).  The  frequency  of  the  pulse-width 
modulation  (PWM)  that  drives  each  boost  converter  is  25  kHz. 

The  solar  cell  converter  (800  W)  is  composed  of  two-phase 
boost  converters  connected  in  parallel  by  the  interleaving  tech¬ 
nique  and  is  selected  to  adapt  the  low  d.c.  voltage  delivered  by  the 
solar  panels,  which  is  approximately  25  V  at  rated  power,  to  the 
60  V  d.c.-bus.  Thus,  each  boost  converter  is  composed  of  a  high- 
frequency  input  inductor  (140  fxH),  four  output  filtering  capacitors 
connected  in  parallel  (total  capacitance:  4400  p,F),  a  diode  (Schot¬ 
tky  diode  RURG3020: 200  V,  30  A),  and  a  main  power  switch  (power 
metal-oxide-semiconductor  field-effect  transistor  IRFP264N  MOS¬ 
FET  :  250  V,  38  A).  The  frequency  of  the  PWM  that  drives  each  boost 
converter  is  25  kHz. 

The  supercapacitor  converter  is  composed  of  four-phase  bidi¬ 
rectional  converters  connected  in  parallel  by  the  interleaving 
technique.  Thus,  each  two-quadrant  converter  is  composed  of 
a  high-frequency  input  inductor  (290  p,H),  an  output  filtering 
capacitor  (aluminum  electrolytic  capacitors:  680  pT),  two  power 
switches  (power  metal-oxide-semiconductor  field-effect  transis¬ 
tor  IRFP264N  MOSFET :  250  V,  38  A),  and  two  diodes  (Schottky  diode 
RURG3020:  200  V,  30  A).  The  frequency  of  the  PWM  that  drives 
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Fig.  5.  Hybrid  power  source  test  bench. 
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Table  1 

d.c.-bus  energy  control  loop  parameters. 


Parameter 

Value 

t'BusREF 

60 

V 

Cbus 

12,200 

pT 

*11 

450 

rad  s-1 

I<n 

22,500 

rad  s~2 

rFC 

0.14 

Q 

rSol 

0.12 

Q 

rsc 

0.10 

Q 

^SCMax 

32 

V 

^SCMin 

15 

V 

fsCRated 

150 

A 

each  bidirectional  converter  is  25  kHz.  The  supercapacitor  current, 
which  flows  across  the  storage  device,  can  either  be  positive  or 
negative,  which  allows  energy  to  be  transferred  in  both  directions. 

3.2.  Control  description 

Measurements  of  the  fuel  cell  current  z'fc,  the  solar  cell  current 
zSol,  the  supercapacitor  current  zSo  the  load  current  zLoad,  the  d.c. 
link  voltage  pbus,  the  fuel  cell  voltage  pFc>  the  solar  cell  voltage  pSoi, 
and  the  supercapacitor  voltage  vSc  are  collected  with  zero-flux  Hall 
effect  sensors. 

The  fuel  cell,  the  solar  cell  and  the  supercapacitor  current- 
regulation  loops  were  realized  by  analog  circuits  to  function  at  a 
high  bandwidth.  Parameters  associated  with  the  d.c.-bus  energy 
regulation  loop  and  the  supercapacitor  energy  regulation  loop  can 
be  seen  in  Tables  1  and  2,  respectively.  The  fuel  cell  power  dynamic 
delay  is  shown  in  Table  2;  this  value  was  experimentally  deter¬ 
mined  as  the  highest  power  slope  of  the  FC  system,  where  no 
fuel  starvation  occurs.  It  must  be  noted  here  that,  for  the  small- 


Table  2 

Supercapacitor  energy  control  loop  parameters. 


Parameter 

Value 

t'SCREF 

25 

V 

Csc 

100 

F 

* 21 

0.1 

wj-1 

PsoiMax  (rated) 

800 

W 

PSolMin 

0 

W 

fsoiMax  (rated) 

28 

A 

JsolMin 

0 

A 

£2 

1 

U)n2 

0.8 

rad  s-1 

PFCMax 

500 

W 

PFCMin 

0 

W 

fFcMax  (rated) 

46 

A 

fFCMin 

0 

A 

1 

COn\ 

0.4 

rads-1 

test  bench,  the  fuel  cell  maximum  power  PFCMax  is  set  at  500  W; 
but,  the  rated  fuel  cell  power  considered  here  is  1200  W.  Further, 
these  two  energy  control  loops,  which  generate  current  refer¬ 
ences  z'fcref.  isoiREF  and  z'scref  are  implemented  on  the  real  time 
card  dSPACE  DS1104  through  the  mathematical  environment  of 
Matlab-Simulink  with  a  sampling  frequency  of  25  kHz. 

3.3.  Experimental  results 

The  studied  d.c.-bus  of  60  V  is  only  connected  to  an  electronic 
load.  To  validate  the  dynamics  of  the  power  (current)  regulation 
loops  for  each  power  source,  Figs.  6-8  present  waveforms  that 
are  obtained  during  the  stepped  power  demand.  Fig.  6  illustrates 
the  fuel  cell  power  demand,  power  response,  voltage  and  current. 
Fig.  7  contains  the  photovoltaic  power  demand,  power  response, 
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Fig.  6.  Dynamic  identification  of  fuel  cell  power  control  loop. 
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Fig.  7.  Dynamic  identification  of  photovoltaic  power  control  loop. 

voltage  and  current.  Finally,  Fig.  8  illustrates  the  supercapacitor 
power  demand,  power  response,  voltage  and  current.  Because  the 
fuel  cell  and  photovoltaic  powers  are  intentionally  limited  by  the 
control  algorithm  (refer  to  Fig.  4:  ‘solar  cell  dynamic  power  limita¬ 
tion’  and  ‘fuel  cell  power  dynamic  limitation’),  a  smooth  transition 
of  the  power  responses  of  the  fuel  cell  and  photovoltaic  sources  is 
observed,  while  the  supercapacitor  power  dynamics  are  very  fast. 
The  supercapacitor  can  supply  power  from  0  to  400  W  in  50  ms.  It 
is  clear  that  the  fast  response  of  the  supercapacitor  storage  device 
can  be  operated  with  the  fuel  cell  and  solar  cell  main  generators 
in  order  to  improve  system  performance.  The  data  in  Fig.  6  also 
confirm  that  the  fuel  cell  is  controlled  to  avoid  the  fuel  starvation 
phenomenon. 

Because  flatness-based  control  is  model-based,  it  may  have 
some  sensitivity  to  errors  in  the  model  parameters.  To  authenti¬ 
cate  its  robustness,  the  flatness-based  control  was  tested  with  the 
exact  model  parameters  (rFC  =  0.14  £2,  rSol  =  0.1 2  £2,  rSc  =  0.1 0  Q)  and 
in  the  lossless  parameters  case  (rFC  =  0  £2,  rSol  =  0  £2,  rSc  =  0  £2).  In  the 
case  of  no  losses,  this  leads  to  a  drastically  less  complex  estima¬ 
tion  of  the  control  law  (Eqs.  (8)— (1 1 )).  Comparisons  (robustness) 
between  the  accurate  parameters  and  the  error  parameters  are 
given  in  Fig.  9.  They  generate  waveforms  that  are  obtained  during 
the  large  load  step  from  0  to  500  W  and  show  the  d.c.-bus  voltage, 
the  load  power  (disturbance),  the  supercapacitor  power,  and  the 
supercapacitor  voltage  (storage  50C).  Because  the  power  of  the  fuel 
cell  and  the  solar  cells  are  intentionally  limited,  the  supercapaci¬ 
tor  supplies  the  transient  load  power  demand.  Similar  waveforms 
are  seen  in  Fig.  9(a)  and  (b).  The  d.c.-bus  voltage  (d.c.  link  stabi¬ 
lization)  is  minimally  influenced  by  the  large  step  in  load  power. 


Clearly  the  performance  of  the  control  system  is  hardly  affected 
by  the  considered  error  in  model  parameters.  Experimental  testing 
demonstrates  that  errors  in  these  parameters  have  relatively  little 
effect  on  regulation  performance.  It  is  therefore  conclude  that  the 
non-linear  differential  flatness-based  approach  provides  an  abso¬ 
lutely  robust  controller  in  this  application. 

Waveforms  obtained  during  the  large  load  cycle  are  presented 
in  Fig.  10.  The  data  show  the  d.c.-bus  voltage,  the  fuel  cell  volt¬ 
age,  the  photovoltaic  voltage,  the  load  power,  the  supercapacitor 
power,  the  fuel  cell  power,  the  photovoltaic  power,  the  superca¬ 
pacitor  current,  the  fuel  cell  current,  the  photovoltaic  current,  and 
the  supercapacitor  voltage  (or  the  supercapacitor  SOC).  In  the  ini¬ 
tial  state,  the  small  load  power  is  equal  to  100  W,  and  the  storage 
device  is  fully  charge,  i.e.,  vsc  =  25  V;  as  a  result,  both  the  fuel  cell 
and  supercapacitor  powers  are  zero,  and  the  photovoltaic  source 
supplies  power  for  the  load  of  100W.  At  t  =  40s,  the  load  power 
steps  to  the  final  constant  power  of  around  900  W  (positive  load 
power  transition).  The  following  observations  are  made: 

•  The  supercapacitor  supplies  most  of  the  900  W  power  that  is 
required  during  the  transient  step  load. 

•  Simultaneously,  the  photovoltaic  power  increases  with  limited 
dynamics  to  a  maximum  power  point  (MPP)  of  around  300  W, 
which  is  limited  by  the  maximum  power  point  tracker  (MPPT). 

•  Concurrently,  the  fuel  cell  power  increases  with  limited  dynamics 
to  a  maximum  power  of  500  W. 

•  The  input  from  the  supercapacitor,  which  supplies  most  of  the 
transient  power  that  is  required  during  the  stepped  load,  slowly 
decreases  and  the  unit  remains  in  a  discharge  state  after  the  load 
step  because  the  steady-state  load  power  (approximately  900  W) 
is  greater  than  the  total  power  supplied  by  the  fuel  cell  and  pho¬ 
tovoltaic  array. 

At  t=100s,  the  supercapacitor  voltage  is  equal  to  19  V.  As  a 
result,  the  supercapacitor  supplies  its  stored  energy  ySc  to  the  d.c.- 
bus.  This  energy  ySCSupply  is  estimated  to  be: 

ysc_suppiy  =  jCSct|c(t  =  40  s)  -  lcsci|c(t  =  100  s)  (27) 

=  13.20  kj 

The  load  power  is  reduced  from  the  high  constant  power  of  900  W 
to  the  low  constant  power  of  100  W  (negative  load  power  transi¬ 
tion).  As  a  result,  the  supercapacitor  changes  from  discharging  to 
charging  and  demonstrates  the  following  four  phases: 

•  First,  the  fuel  cell  and  photovoltaic  array  still  supply  their  con¬ 
stant  maximum  powers  to  drive  the  load  and  to  charge  the 
supercapacitor. 

•  Second,  at  t=  1 10  s,  the  supercapacitor  is  approaching  full  charge, 
i.e.,  vSc  =  23  V.  Consequently,  the  fuel  cell  power  is  reduced  with 
limited  power  dynamics. 

•  Third,  at  t=  120  s,  the  supercapacitor  is  nearly  fully  charged,  i.e., 
vsc  =  24  V.  As  a  result,  the  photovoltaic  power  is  reduced  with 
limited  power  dynamics. 

•  Fourth,  at  t  =  160  s,  the  supercapacitor  is  fully  charged,  i.e.,  vSc  = 
^scref  =  25  V.  After  slowly  decreasing,  the  photovoltaic  power 
remains  at  a  constant  power  of  100W  for  the  load  power 
demanded.  Furthermore,  the  fuel  cell  and  supercapacitor  powers 
are  zero. 

It  is  evident  that  the  d.c.-bus  voltage  waveform  is  stable  during 
the  large  load  cycle,  which  is  critically  important  when  employing 
supercapacitors  to  improve  the  dynamic  performance  of  the  whole 
system  using  the  proposed  control  law. 

Finally,  Fig.  1 1  presents  waveforms  that  are  obtained  during  the 
long  load  cycles  measured  on  5  December,  2009.  The  waveforms  are 
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Fig.  8.  Dynamic  identification  of  supercapacitor  power  control  loop. 
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Fig.  9.  Comparison  of  d.c.  link  stabilization  of  hybrid  energy  source  during  large  load  step,  (a)  exact  model  (rFc  =  0.14^2,  rSol  =  0.12£2,  rSc  =  0.10^2)  and  (b)  error  model 
(robustness)  (rFc  =  0  £2,  rSoi  =  0  £2,  rSc  =  0  £2). 
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Fig.  10.  Hybrid  source  response  during  load  cycle. 


esults  in  Fig.  10.  During  the  experiment, 
i/vas  set  to  500  W,  and  the  PV  maximum 
r  radiation,  weather  conditions  and  tem- 
ae  MPPT.  For  example,  at  1 5 : 01 :40,  the  PV 
iroximately  equal  to  180  W;  at  15:02:30, 
ras  about  1 80  W;  at  1 5 :03 :20,  the  PV max- 
180  W;  and  from  15:05:50  to  15:08:20, 
/as  reduced  to  0  W,  because  of  the  cloudy 


conditions.  In  particular,  it  was  found  that  the  power  plant  was 
always  energy  balanced  pLoad(£) » PfcM + PsoiW + Psc(f)  by  the  pro¬ 
posed  original  control  algorithm. 

The  important  variable  necessary  to  balance  the  energy  in  this 
complex  system  is  the  d.c.-bus  energy  or  voltage.  From  the  exper¬ 
imental  validation,  the  d.c.-bus  voltage  is  automatically  controlled 
at  the  constant  set-point,  i.e.,  pbus  =  60  V.  This  experiment  confirms 
that  the  energy  in  the  system  is  well  managed.  The  fuel  cell  and 
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Fig.  11.  Hybrid  source  response  during  long  load  cycles. 


solar  cell  powers  are  limited  at  their  maximum  powers,  and  the 
fuel  cell  power  dynamics  are  controlled;  as  a  result,  there  is  no 
fuel  starvation  problem  and  the  fuel  cell  stack  lifetime  is  increased 
[46-48]. 

4.  Conclusions 

Energy  management  of  multi-power  sources  has  been  proposed 
as  a  solution  for  a  hybrid  energy  system  that  uses  renewable  energy 
from  solar  cells,  fuel  cells  and  a  supercapacitor  as  an  energy  stor¬ 


age  device.  A  supercapacitor  can  advance  the  load,  following  the 
characteristics  of  the  main  sources  by  providing  a  stronger  power 
response  to  changes  in  the  system  load.  During  essential  steps  in  the 
load,  the  supercapacitor  provides  the  energy  balance  needed  dur¬ 
ing  load  transition  periods.  Adding  energy  storage  to  the  distributed 
power  systems  improves  power  quality  and  efficiency. 

Experimental  verification  with  a  small-scale  hybrid  test  bench 
(Nexa  Ballard  fuel  cell  power  generator:  1.2  kW,  46  A;  Ekarat  Solar 
Cell  power  module:  800  W,  31  A;  Maxwell  supercapacitor  storage 
device:  1 00  F,  32  V)  has  demonstrated  the  excellent  performance  of 
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the  whole  system,  and  has  validated  the  proposed  energy  manage¬ 
ment  principle. 
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